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INTRODUCTION 


X 


Since  the  earliest  use  of  radar,  radar  operators  have  detected  and  tracked  targets  using  visual 
inputs  from  PPIs  and  A-scopes.  Although  operators  can  perform  these  tasks  very  accurately,  they  are 
easily  saturated  and  quickly  become  fatigued.  To  correct  this  situation,  automatic  detection  and  track¬ 
ing  (ADT)  systems  are  being  attached  to  many  radars.  Undoubtedly,  as  digital  processing  increases  in 
speed  and  hardware  decreases  in  cost  and  size,  ADT  systems  will  become  associated  with  almost  all  but 
the  simplest  radars. 


The  purpose  of  this  report  is  to  survey  the  state  of  the  art  in  ADT  and  update  the  ADT  portion 
included  in  "Survey  of  Radar  Signal  Processing"  [11.  Included  in  this  discussion  are  various  non¬ 
coherent  integrators  that  provide  target  enhancement,  thresholding  techniques  for  reducing  false  alarms 
and  target  suppression,  and  algorithms  for  estimating  target  position  and  resolving  targets.  The  other 
major  topic  that  is  covered  is  track-while-scan  systems.  An  overview  of  the  entire  tracking  system  is 
given,  followed  by  a  discussion  of  its  various  components  such  as  tracking  filter,  maneuver-following 
logic,  track  initiation,  and  correlation  logic.  Finally,  the  survey  concludes  with  a  discussion  of  radar 
netting,  including  both  colocated  and  multisite  systems. 

NONCOHERENT  INTEGRATORS 


There  are  many  different  integrators  that  are  used  with  radars.  While  the  optimal  integrator  is 
never  implemented,  it  will  be  introduced  in  order  to  compare  the  various  practical  integrators  to  it. 

Optimal  Detector 


The  radar  detection  problem  is  a  binary  hypothesis  testing  problem;  i.e., 

H0-  No  target  present 

//j:  Target  present 

While  many  criteria  can  be  used  to  solve  this  problem,  the  most  appropriate  for  radar  is  the  Nsyman- 
Pearson  criterion.  This  criterion  maximizes  the  probability  of  detection  PD  for  a  given  probability  of 
false  alarm  Pfa  by  comparing  the  likelihood  ratio  to  an  appropriate  threshold.  It  can  be  shown  that  the 
likelihood  ratio  can  be  reduced  so  that  one  should  declare  a  target  present  in  the  j th  range  cell  if 


n  A-X- 

7 r  70  >  r. 

i-i  <r 


where  n  is  the  number  of  pulses  within  the  beamwidth,  xu  is  the  /th  envelop-detected  returned  pulse  in 
the  ,/th  range  cell.  A,  is  the  signal  amplitude  on  the  /th  pulse  ({/4,}  are  proportional  to  the  square  of  the 
voltage  antenna  pattern  as  the  antenna  sweeps  over  the  target),  a-2  is  the  noise  power,  70  is  the  Bessel 
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function  of  zero  order,  and  the  threshold  T  determines  the  Pfa-  For  small  signals,  the  optimal  detector 
can  be  approximated  by 

t  *}  4  >  r. 

i~l 

Most  of  the  detectors  implemented  differ  from  the  optimal  detector  in  the  weights  (e.g..  A,)  applied  to 
the  pulses. 

Besides  detecting  targets,  the  detection  process  involves  making  angular  estimates  of  the  azimuth 
position  of  the  target.  Swerling  12]  calculated  the  standard  deviation  of  the  optimal  estimate  by  using 
the  Cramer-Rao  bound.  His  result  holds  for  a  large  number  of  pulses  within  the  beamwidth  and  the 
optimal  estimate  involves  finding  the  location  where  the  correlation  of  the  returned  signal  and  the 
derivative  of  the  fourth  power  of  the  voltage  antenna  pattern  is  zero.  Though  this  estimate  is  rarely 
implemented,  its  performance  is  approached  by  simple  estimates  such  as  the  maximum  value  of  the 
integrator  or  the  midpoint  between  the  initial  and  final  crossing  of  the  detection  threshold. 

A  few  of  the  most  common  integrators  are  shown  in  Fig.  1.  Though  they  are  shown  as  being 
implemented  with  shift  registers,  with  today’s  technology  they  would  probably  be  implemented  with 
random  access  memory.  The  input  to  these  detectors  can  be  either  linear,  squ?re-law,  or  log  video. 
Since  linear  video  is  probably  the  most  commonly  used  video,  the  advantages  ana  disadvantages  of  the 
various  integrators  will  be  stated  for  this  video. 
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Moving  Window 

The  moving  window  performs  a  running  sum  of  n  pulses  in  each  range  cell; 

S,  -  S,_!  +  x,  - 

where  S,  is  the  sum  on  the  /th  pulse  and  for  notational  convenience  the  subscript  for  the  y'th  range  cell 
has  been  dropped.  The  detection  performance  of  this  detector  is  only  0.5  dB  worse  than  the  optimal 
detector.  The  angular  estimate  that  is  obtained  by  either  taking  the  maximum  value  of  the  running 
sum  or  taking  the  midpoint  between  the  first  and  last  crossing  of  the  detection  threshold  has  a  bias  of 
n/2  pulses  which  is  easily  corrected.  The  standard  deviation  of  the  estimation  error  of  both  estimators 
is  about  20%  higher  than  the  Cramer-Rao  lower  bound.  The  major  disadvantage  of  this  detector  is  that 
the  last  n  pulses  for  each  range  cell  must  be  saved,  resulting  in  a  large  storage  requirement  for  large 
beamwidth  radars.  However,  because  of  the  current  availability  of  large  memories  of  reduced  size  and 
cost,  this  problem  has  essentially  disappeared. 

Feedback  Integrator 

The  amount  of  storage  required  can  be  reduced  significantly  by  using  a  feedback  integrator 

S-KS^+x,. 

Although  the  feedback  integrator  applies  an  exponential  weighting  into  the  past,  its  detection  perfor¬ 
mance  is  only  1  dB  less  than  the  optimal  integrator.  Unfortunately,  difficulties  are  encountered  when 
using  the  feedback  integrator  to  estimate  the  azimuth  position.  The  threshold  crossing  procedure  yields 
estimates  only  20%  greater  than  the  lower  bound,  but  the  bias  is  a  function  of  the  signal-to-noise  ratio 
(S/N)  and  must  be  estimated.  On  the  other  hand,  the  maximum  value,  though  it  has  a  constant  bias, 
has  estimates  that  are  100%  greater  than  the  lower  bound.  Furthermore,  the  exponential  weighting 
function  essentially  destroys  the  radar  sidelobes.  Because  of  these  problems,  I  would  recommend  that 
the  feedback  integrator  never  be  implemented. 

Two-Pole  Filter 

The  two-pole  filter  requires  the  storage  ol  an  intermediate  calculation  in  addition  to  the  integrated 
output  (see  Fig.  1).  However,  with  this  rather  simple  device,  a  weighting  pattern  similar  to  the  antenna 
pattern  can  be  obtained.  The  detection  performance  is  within  0.15  dB  of  the  optimal  detector  (i.e., 
better  than  the  moving  window)  and  its  angular  estimates  are  about  20%  greater  than  the  Cramer-Rao 
lower  bound.  The  problems  with  this  detector  are  1.  it  has  rather  high  detector  sidelobes,  15  to  20  dB 
and  2,  it  is  extremely  sensitive  to  interference  (i.e.,  the  filter  has  a  high  gain  resulting  in  a  large  output 
for  a  single  sample  that  has  a  high  value). 

Binary  Integrator 

The  binary  integrator  is  also  called  a  binary  moving  window,  an  w-out-of-n  detector,  a  dual  thres¬ 
hold  detector,  or  rank  detector.  The  input  samples  are  quantized  to  0  or  1  depending  on  whether  or 
not  they  are  greater  than  a  threshold  T\.  The  last  n  0s  and  Is  ire  summed  and  compared  to  a  detection 
threshold  F2.  For  a  large  number  of  pulses,  the  detection  performance  of  this  detector  is  2  dB  less  than 
the  moving  window  because  of  the  hard  limiting  of  the  data,  and  the  angular  estimation  error  is  25% 
greater  than  the  Cramer-Rao  lower  bound.  There  are  several  reasons  why  this  detector  is  used:  1.  it  is 
easily  implemented,  2.  it  ignores  interference  spikes  which  cause  trouble  with  integrators  that  directly 
use  signal  amplitude,  and  3.,  it  works  extremely  well  when  'he  noise  has  a  non-Rayleigh  density.  For 
log-normal  noise  statistics,  the  binary  integrator  is  less  than  a  decibel  from  the  optimal  detector  while 
the  moving  window  is  10  dB  from  the  optimal  detector  [3].  The  binary  integrator  should  always  be 
used  if  there  are  only  two  pulses  on  target,  a  condition  that  occurs  with  many  3-D  radars. 
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Batch  Processor 

This  detector  is  a  combination  of  an  integrator  and  dump  and  a  binary  integrator  (it  could  also  use 
a  moving  window  instead  of  a  binary  integrator).  The  detector  should  be  used  when  there  are  a  large 
number  of  pulses.  For  instance,  if  there  are  80  pulses  within  the  3-dB  beamwidth,  one  could  integrate 
16  pulses  (this  is  called  a  batch)  and  then  use  the  five  batches  as  an  input  to  either  a  binary  integrator 
or  moving  window  integrator.  The  performance  is  only  slightly  worse  than  the  moving  window  but  the 
batch  integrator  is  much  easier  to  implement  than  the  moving  window.  The  angle  estimate  is  given  by 

-  £  Bfl, 

9°  £fl,  ' 

where  B,  is  batch  amplitude  and  9,  is  azimuth  angle  corresponding  to  the  center  of  the  batch. 

FALSE  ALARM  CONTROL 

If  fixed  thresholds  are  used  with  the  previously  discussed  integrators  the  detections  will  resemble 
the  output  of  a  PPI.  In  the  presence  of  clutter,  this  can  result  in  an  enormous  number  of  detections 
which  will  saturate  and  disrupt  the  tracking  computer  associated  with  the  radar  system.  Of  course,  one 
can  limit  the  number  of  false  alarms  with  a  fixed  threshold  system  by  setting  a  very  high  threshold. 
Unfortunately,  this  would  reduce  target  sensitivity  in  regions  of  low  noise  (clutter)  return.  What  is 
desired  is  a  detector  that  will  detect  a  target  when  it  has  a  higher  return  than  its  immediate  background. 
Two  such  approaches  are  adaptive  thresholding  and  nonparametric  detectors.  Both  of  these  approaches 
assume  that  the  samples  in  the  range  cells  surrounding  the  test  cell  (called  reference  cells)  are  indepen¬ 
dent  and  identically  distributed;  and  furthermore,  it  is  usually  assumed  that  the  time  samples  are 
independent.  Both  kinds  of  detectors  test  whether  the  test  cell  has  a  return  sufficiently  larger  than  the 
reference  cells. 

Adaptive  Thresholding 

The  basic  assumption  of  the  adaptive  thresholding  technique  is  that  the  noise  density  is  known 
except  for  a  few  unknown  parameters.  The  surrounding  reference  cells  are  then  used  to  estimate  the 
unknown  parameters  and  a  threshold  based  on  the  estimated  parameters  is  obtained.  The  simplest 
adaptive  detector,  shown  in  Fig.  2,  is  the  cell-averaging  CFAR  (constant  false  alarm  rate)  investigated 
by  Finn  [4J.  If  the  noise  has  a  Rayleigh  density,  only  the  parameter  o-  needs  to  be  estimated  and  the 
threshold  is  of  the  form  T  «  Klx,  —  Knyfirjl  a.  However,  since  T is  set  by  an  estimate  <r,  it  must  be 
slightly  larger  than  the  threshold  one  would  use  if  cr  were  known  a  priori.  The  raised  threshold  causes 
a  loss  in  target  sensitivity  and  is  referred  to  as  a  CFAR  loss.  This  loss  has  been  calculated  and  is  given 
in  (5]. 

If  there  is  uncertainty  about  whether  or  not  the  noise  is  Rayleigh  distributed,  it  is  better  to  thres¬ 
hold  the  data  earlier  and  use  a  binary  integrator  as  shown  in  Fig.  3.  The  reason  why  this  detector  is 
robust  to  the  noise  density  is  that  by  setting  K  to  yield  a  "1"  with  probability  0.1,  a  Pfa  =  10~6  can  be 
obtained  by  using  a  7-out-of-9  detector.  While  noise  may  be  non-Rayleigh,  it  will  probably  be  very 
Raleigh-like  out  to  the  tenth  percentile.  Furthermore,  one  can  use  a  slow  feedback  to  control  K  in 
order  to  maintain  a  desired  Pfa  either  on  a  scan  or  sector  basis. 

If  the  noise  power  varies  from  pulse-to-pulse  basis  (as  it  would  in  jamming  when  frequency  agility 
is  employed)  one  must  CFAR  each  pulse  and  then  integrate.  While  the  binary  integrator  performs  this 
type  of  CFAR  action,  analysis  [6]  has  shown  that  the  ratio  detector  shown  in  Fig.  4  is  a  better  detector. 
A  comparison  of  the  ratio  detector  with  other  commonly  used  detectors  when  the  jamming  level  varies 
by  20  dB  per  pulse  is  shown  in  Fig.  5. 
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Fig.  1  —  Cell-averaging  CFAR 
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Fig.  3  —  l  -"plementation  of  binary  integrator 


Fig.  4  —  Ratio  detector 


NRL  REPORT  8698 


Clutter  Mapping 

A  clutter  map  uses  adaptive  thresholding  where  the  threshold  is  calculated  from  the  return  in  the 
test  cell  on  previous  scans  rather  than  from  the  surrounding  reference  cells  on  the  same  scan.  This 
technique  has  the  advantage  that  for  essentially  stationary  environments  (e.g.,  land  based  radar  against 
ground  clutter),  the  radar  has  inter  clutter  visibility— it  can  see  between  large  clutter  returns.  Lincoln 
Laboratory  [7]  in  their  Moving  Target  Detector  (MTD)  used  a  clutter  map  (setting  their  threshold 
using  the  return  from  the  last  eight  scans)  for  the  zero  doppler  filter  very  effectively.  There  remain 
several  unanswered  questions  concerning  the  use  of  clutter  maps  with  radars  whose  frequency  agility 
extends  over  a  large  bandwidth. 
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Nonparametric  Detectors 

The  most  common  way  nonparametric  detectors  obtain  CFAR  is  by  ranking  the  test  sample  with 
the  reference  cells.  Under  the  hypothesis  that  all  the  samples  are  independent  samples  from  an  un¬ 
known  density  function,  the  test  sample  has  a  uniform  density  function.  For  instance,  referring  to  the 
ranker  in  Fig.  6,  the  test  cell  is  compared  to  15  of  its  neighbors.  Since  in  the  set  of  16  samples,  the  test 
sample  has  equal  probability  of  being  the  smallest  sample  (or  equivalently  any  other  rank),  the  proba¬ 
bility  that  the  test  sample  takes  on  values  0,  1,  ....  15  is  1/16.  A  simple  rank  detector  is  constructed 
by  comparing  the  rank  to  a  threshold  K\  and  putting  out  a  1  if  the  rank  is  larger,  a  0  otherwise.  The  Os 
and  Is  are  summed  in  a  moving  window.  This  detector  incurs  a  CFAR  loss  of  about  2  dB,  but  achieves 
a  fixed  Pfa  for  any  unknown  noise  density  as  long  as  the  time  samples  are  independent.  This  detector 
is  incorporated  into  ARTS-3A  post  processor  used  in  conjunction  with  the  FAA’s  ASR  radar.  The 
shortcomings  of  this  detector  are  1.  it  is  fairly  susceptible  to  target  suppression  (e.g.,  if  a  large  target  is 
in  the  reference  cells,  the  test  cell  cannot  receive  the  highest  ranks)  and  2.  only  certain  values  of  Pfa 
can  be  obtained. 


K 

Fig.  6  —  Rant  detector 


If  the  time  samples  are  dependent,  the  rank  detector  will  not  yield  CFAR.  A  modified  rank 
detector,  called  the  modified  generalized  sign  test  (MGST)  [8),  is  an  attempt  to  maintain  a  low  Pfa  and 
is  shown  in  Fig.  7.  This  detector  can  be  divided  into  three  parts:  a  ranker,  an  integrator  (in  this  case  a 
two-pole  filter),  and  a  threshold  (decision  process),  A  target  is  declared  when  the  integrated  output 
exceeds  two  thresholds.  The  first  threshold  is  fixed  (equals  n  +  T\/K  from  Fig.  7)  and  yields 
Pfa  =  10  6  when  the  reference  cells  are  independent  and  identically  distributed.  The  second  threshold 
is  adaptive  and  maintains  a  low  Pfa  when  the  reference  samples  are  correlated.  The  device  estimates 
the  standard  deviation  of  the  correlated  samples  with  the  mean  deviate  estimator,  where  extraneous  tar¬ 
gets  in  the  reference  cells  have  been  excluded  from  the  estimate  by  use  of  a  preliminary  threshold  r2. 

The  rank  and  MGST  detectors  are  basically  2-sample  detectors.  They  decide  a  target  is  present  if 
the  ranks  of  the  test  cell  are  significantly  greater  than  the  ranks  of  the  reference  cells.  Target  suppres¬ 
sion  occurs  at  all  interfaces  (e.g.,  land,  sea),  where  the  homogeneity  assumption  is  violated.  However, 
there  exist  some  tests,  such  as  Spearman  Rho  and  Kendall  Tau  tests,  that  only  depend  on  the  test  cell. 
These  tests  work  on  the  fact  that  as  the  antenna  beam  sweeps  over  a  point  target,  the  signal  return 
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increases  and  then  decreases.  Thus,  for  the  test  cell,  the  ranks  should  follows  a  pattern— first  increas¬ 
ing  and  then  decreasing.  Although  these  detectors  do  not  require  reference  cells  and  hence  have  the 
useful  nroperty  of  not  requiring  homogeneity,  these  detectors  are  not  generally  used  because  of  the 
large  CFAR  loss  taken  for  moderate  sample  sizes  [9]. 


A  basic  disadvantage  of  all  nonparametric  detectors  is  that  one  loses  amplitude  information,  which 
can  be  a  very  important  discriminant  between  target  and  clutter  [10].  For  example,  a  large  return  ( <r  > 
1000  m2)  in  a  clutter  area  is  probably  just  clutter  breakthrough. 


Target  Suppression 


All  detectors  are  susceptible  to  target  suppression  to  one  extent  or  another.  There  are  basically 
three  approaches  to  this  problem:  one  can  use  clutter  maps  which  have  no  reference  cells;  secondly, 
one  can  try  to  remove  large  returns  from  the  calculation  of  the  threshold  [8, 11];  or  finally,  one  can  try 
to  diminish  the  effects  of  large  returns  by  either  limiting  or  using  log  video.  The  technique  that  should 
be  used  is  a  function  of  the  particular  radar  system  and  its  environment. 


Target  Resolution 


A  single,  large  target  will  probably  be  detected  many  times,  e.g.,  in  adjacent  range  cells,  azimuth 
beams,  and  elevation  beams.  Therefore,  all  automatic  detectors  have  algorithms  for  merging  the  indivi¬ 
dual  detections  into  a  centroided  detection.  However,  most  algorithms  have  been  designed  so  that  they 
will  rarely  ever  split  a  single  target  into  two  targets— a  decision  resulting  in  poor  range  resolution  capa¬ 
bility.  For  instance,  a  commonly  used  merging  algorithm  [12]  has  a  range  resolution  of  approximately 
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two  and  a  half  pulsewidths.  Inherently,  the  range  resolution  is  a  pulsewidth  or  less  depending  on  the 
signal-to-noise  ratios  and  the  phase  difference  between  the  two  targets.  Since  the  number  of  targets  in 
a  raid  can  be  important,  this  problem  is  currently  under  investigation. 

TRACKING  SYSTEM 

A  general  outline  of  a  track-while-scan  system  (tracking  systems  for  surveillance  radars  whose 
nominal  scan  time  is  from  4  to  12  s)  is  considered  first.  Then  the  tracking  filter,  maneuver-following 
logic,  tracking  initiation,  and  correlation  logic  is  discussed  in  detail.  Finally,  methods  of  integrating  data 
from  several  radars  is  discussed. 

System  Outline 

Almost  all  track-while-scan  systems  operate  on  a  sector  basis.  A  typical  series  of  operations  is 
shown  in  Fig.  8.  For  instance,  if  the  radar  has  reported  all  the  detections  in  sector  11  and  is  now  in 
sector  12,  the  tracking  program  would  start  by  correlating  (trying  to  associate)  the  clutter  points  (sta¬ 
tionary  tracks)  in  sector  10  with  detections  in  sectors  9,  10,  and  11.  Those  detections  that  are  associ¬ 
ated  with  clutter  points  are  deleted  (are  not  used  for  further  correlations)  from  the  detection  file  and 
are  used  to  update  the  clutter  points.  Next,  firm  tracks  in  sector  8  are  correlated  with  detections  in  sec¬ 
tors  7.  8,  and  9.  By  this  time  all  clutter  points  have  been  removed  from  sectors  9  and  below.  Those 
detections  that  are  associated  with  firm  tracks  are  deleted  from  the  detection  file  and  are  used  to  update 
the  appropriate  track. 


Fig.  8  —  Various  operations  of  a  track-white-scan  system 
performed  on  a  sector  basis 


Usually,  some  provision  is  made  for  giving  preference  to  firm  tracks  (with  respect  to  tentative 
tracks)  in  the  correlation  process.  By  performing  the  correlation  process  two  sectors  behind  firm  track 
correlations  (Fig.  8),  it  is  impossible  for  tentative  tracks  to  steal  detections  belonging  to  firm  tracks.  In 
other  tracking  systems  the  correlation  for  firm  and  tentative  tracks  is  performed  in  the  same  sector, 
however,  the  generalized  distance  D  between  tracks  and  detections  is  incremented  by  A D  if  the  track  is 
tentative. 

Finally,  detections  that  are  not  associated  with  either  clutter  points  or  tracks  are  used  for  initiation 
purposes.  The  most  common  initiation  procedure  is  to  initiate  a  tentative  track;  later  the  track  is 
dropped  or  else  made  a  firm  track  or  clutter  point.  An  alternate  approach  [13]  is  to  establish  both  a 
clutter  point  and  a  tentative  track.  If  the  detection  came  from  a  stationary  target,  the  clutter  point  will 
be  updated  and  the  tentative  track  will  eventually  be  dropped.  On  the  other  hand,  if  the  detection  came 
from  a  moving  target,  the  tentative  track  will  be  made  firm  and  the  clutter  point  will  be  dropped.  The 
latter  method  requires  less  computer  computation  time  when  most  of  the  detections  are  clutter  resi¬ 
dues. 
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The  correlation  procedure  is  made  in  a  sector  framework  to  avoid  the  necessity  of  correlating  all 
tracks  with  all  detections.  The  procedure  can  be  implemented  very  easily  by  defining  two  computer 
arrays:  a  sector  file  and  a  track  file.  The  sector  file  for  sector  I  contains  the  first  track  number  in  sector 
I,  and  the  track  file  for  track  /  contains  the  next  track  number  in  the  same  sector  as  track  /  or  a  0,  indi¬ 
cating  that  the  track  is  the  last  track  in  the  sector. 

Tracking  Filters 

Before  proceeding,  the  coordinate  system  in  which  the  tracking  will  be  performed  is  discussed. 
The  target  location  measured  by  the  radar  is  in  spherical  coordinates:  range,  azimuth,  elevation,  and 
possibly  range  rate.  Thus,  it  may  seem  natural  to  perform  tracking  in  spherical  coordinates.  However, 
this  causes  difficulties  since  the  motion  of  constant-velocity  targets  (straight  lines)  will  cause  accelera¬ 
tion  terms  in  all  coordinates.  A  simple  solution  to  this  problem  is  to  track  in  a  cartesian  coordinate  sys¬ 
tem.  Although  it  may  appear  that  tracking  in  cartesian  coordinates  will  destroy  the  accurate  range  track, 
it  has  been  shown  [14]  that  the  inherent  accuracy  is  maintained. 

The  most  commonly  used  filter  is  the  a-0  filter  described  by 

Xjik)  -  x„(k)  +  aUm(fc)  -  ^,(/c)], 

V,(k)  -  Vs(k  -  1)  +  filxjk)  -  xp(k)VT, 
and 

Xp(k  +  1)  —  xs(k)  +  Vs(k)T. 

where  x,(k)  is  the  smoothed  position,  V,(k)  is  the  smoothed  velocity,  Xp(k)  is  the  predicted  position, 
xm(k)  is  the  measured  position,  T  is  the  scanning  period  (time  between  detections),  and  a  and  0  are 
the  system  gains. 

The  optimal  filter  for  performing  the  'racking  when  the  equation  of  motion  is  known  is  the  Kal¬ 
man  filter  [15],  The  Kalman  filter  is  a  recursive  filter  that  minimizes  the  least-square  error.  The  Kal¬ 
man  filter  is  usually  not  implemented  because  of  the  complexity  of  the  calculations  required.  It  is  fairly 
simple  to  show  that  under  certain  conditions  the  a-0  filter  can  be  made  equivalent  to  the  Kalman  filter 
by  setting 

2(2 k  -  1) 

*(*  +  1)  ' 
and 

k(k  +  1) 

on  the  kth  scan.  Thus,  as  time  passes,  a  and  0  approach  zero,  applying  heavy  smoothing  to  the  new 
samples.  In  practice,  a  and  0  usually  take  on  only  three  values:  an  initial  value,  an  intermediate  value, 
and  a  steady-state  value. 

Maneuver-Following  Logic 

Benedict  and  Bordner  [16]  noted  that  in  track-while-scan  systems  there  is  a  conflicting  require¬ 
ment  between  good  noise  smoothing  (implying  small  a  and  0)  and  good  maneuver-following  capability 
(implying  large  a  and  0).  Although  some  compromise  is  always  required,  the  smoothing  equations 
should  be  constructed  to  «ive  the  best  compromise  for  a  desired  noise  reduction.  Benedict  and  Bordner 
defined  a  measure  of  transient-following  capability  and  showed  that  a  and  0  should  be  related  by 
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Thus,  an  (a,/3)  pair  satisfying  the  above  equation  can  be  chosen  so  that  the  tracking  filter  will  follow  a 
specified  g  turn. 

The  trouble  with  the  preceding  method  is  that  if  high-#  turns  are  followed,  the  noise  performance 
is  rather  poor.  To  rectify  this  situation,  a  turn  detector  employing  the  two  correlation  regions  shown  in 
Fig.  9  is  used.  If  the  detection  is  in  the  nonmaneuvering  correlation  region,  the  filter  operates  as  usual. 
When  the  target  falls  outside  the  inner  gate  but  within  the  maneuver  gate,  a  maneuver  is  declared  and 
the  filter  bandwidth  is  increased  (a  and  /i  ai^  increased).  To  avoid  the  problem  of  the  target  fading  and 
a  false  alarm  appearing  in  the  large  maneuver  gate,  the  track  should  be  bifurcated  when  a  maneuver  is 
declared.  That  is,  two  tracks  are  generated:  the  old  track  with  no  detection  and  a  new  maneuvering 
track  with  the  new  detection  and  increased  bandwidth.  The  next  detection  (or  two)  is  used  to  resolve 
the  ambiguity  and  remove  one  of  the  tracks.  Other  solutions  are  to  adjust  the  bandwidth  as  a  function 
of  the  measurement  error  [13]  and  to  use  the  Kalman  filter  with  a  realistic  target-maneuvering  model 
[17]. 


MANEUVERING  GATE 


.  NONMANEUVERING 
GATE 


POSITION 


Fig.  9  —  Maneuver  and  nonmaneuver  gates  centered 
at  the  target’s  predicted  position 


Track  Initiation 

Detections  that  do  not  correlate  with  clutter  points  or  update  tracks  are  used  to  initiate  new  tracks. 
If  the  detection  does  not  contain  doppler  information,  the  new  detection  is  usually  used  as  the  predicted 
position  (in  some  systems,  one  assumes  a  radially  inbound  velocity),  and  a  large  correlation  region 
must  be  used.  Since  the  probability  of  obtaining  false  alarms  in  the  large  correlation  region  is  large, 
one  should  drop  the  track  if  no  correlation  is  obtained  on  the  second  scan;  and  also  one  should  never 
declare  a  track  firm  until  a  third  detection  (falling  within  a  smaller  correlation  region)  is  obtained.  The 
usual  initiation  criteria  are  three  out  of  four  and  three  out  of  five  although  one  may  require  four  detec¬ 
tions  in  regions  with  high  false  alarm  rates.  The  possible  exceptions  for  using  only  two  detections  are 
when  doppler  information  is  available  (so  that  a  small  correlation  region  can  be  used  immediately)  or 
for  "pop-up"  (close)  targets  in  a  military  situation. 

An  alternative  track  initiation  logic  [18]  uses  a  sequential  hypothesis-testing  scheme.  When  a 
correlation  is  made  on  the  ith  scan.  A,  is  added  to  the  likelihood;  when  a  correlation  opportunity  is 
missed,  A/  is  subtracted  from  the  likelihood.  The  increments  are  set  by  the  state  of  the  tracking  system 
and  are  a  function  of  the  closeness  of  the  association,  the  number  of  false  alarms,  the  a  priori  probabil¬ 
ity  of  targets,  and  the  probability  of  detection.  Although  this  method  will  inhibit  false  tracks  in  dense 
detection  environments,  it  will  not  necessarily  establish  the  correct  tracks.  To  initiate  tracks  in  a  dense 
detection  environment,  APL  [19]  has  investigated  "retrospective  processing"  where  the  detections  over 
the  last  eight  scans  are  saved  and  used  to  initiate  straight-line  tracks  using  a  collection  of  filters  matched 
to  different  velocities. 

Firm  tracks  that  are  not  updated  in  30  to  40  seconds  are  usually  dropped. 
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Correlation  Logic 

Several  procedures  will  now  be  given  for  associating  detections  with  tracks.  Of  special  interest  are 
the  conflicting  situations  of  multiple  tracks  competing  for  a  single  detection  or  of  multiple  detections 
lying  within  a  track’s  correlation  gate. 

To  limit  the  number  of  detections  that  can  update  a  track,  correlation  gates  are  used.  A  detection 
can  never  update  a  track  unless  it  lies  within  the  correlation  gate  that  is  centered  at  the  track’s  predicted 
position.  The  correlation  gate  should  be  defined  in  r-0  coordinates,  regardless  of  what  coordinate  sys¬ 
tem  is  being  used  for  tracking.  Furthermore,  the  gate  size  should  be  a  function  of  the  measurement 
accuracy  and  prediction  error  so  that  the  probability  of  the  correct  detection  lying  within  the  gate  is  high 
(at  least  0.99),  In  some  tracking  systems,  the  correlation  gate  is  fed  back  to  the  automatic  detector, 
and  the  detection  threshold  is  lowered  in  the  gate  to  increase  the  PD. 

When  several  detections  are  within  the  correlation  region,  the  usual  and  simplest  solution  is  to 
associate  the  closest  detection  with  the  track.  Specifically,  the  measure  of  closeness  is  the  statistical  dis¬ 
tance 

D2  „  (rp  -  r")2  +  (Op  -  O2 

ffr  <T« 

v'here  ( r 9P)  is  the  predicted  position  (rm,  0m)  is  the  measured  position,  or}  is  the  variance  of 
rP  ~  rnr>  and  o-g  is  the  variance  of  0p  -  0m,  Since  the  prediction  variance  is  proportional  to  the  meas¬ 
urement  variance,  o-g  and  6}  are  sometimes  replaced  by  the  measurement  variances.  Statistical  distance 
rather  than  Euclidean  distance  must  be  used  because  the  range  accuracy  is  usually  much  better  than  the 
azimuth  accuracy. 

Problems  associated  with  multiple  detections  and  tracks  are  illustrated  in  Fig.  10;  two  detections 
are  within  gate  1,  three  detections  are  within  gate  2,  and  one  detection  is  within  gate  3.  Table  1  lists  all 
detections  within  the  tracking  gate,  and  the  detections  are  entered  in  the  order  of  their  statistical  dis¬ 
tance  from  the  track.  Tentatively,  the  closest  detection  is  associated  with  each  track,  and  then  the  ten¬ 
tative  associations  are  examined  to  remove  detections  that  are  used  more  than  once.  Detection  8, 
which  is  associated  with  tracks  1  and  2,  is  paired  with  the  closest  track  (track  1  in  this  case);  then  all 
other  tracks  are  reexamined  to  eliminate  all  associations  with  detection  8.  Detection  7  is  associated 
with  tracks  2  and  3;  the  conflict  is  noted  but  is  resolved  by  pairing  detection  7  with  track  2.  When  the 
other  association  with  detection  7  is  eliminated,  track  3  has  no  associations  with  it  and  consequently  will 
not  be  updated  on  this  scan.  Thus  track  1  is  updated  by  detection  8,  track  2  is  updated  by  detection  7, 
and  track  3  is  not  updated. 
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Fig.  10  —  Examples  of  the  problems  caused  by  multiple 
detections  and  tracks  in  close  vicinity 
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Table  1  —  Association  Table  for  the  Example  Shown  in  Fig.  10 


Track 

Closest  Association 

Second  Association 

Third  Association 

Number 

Detection 

Number 

D2 

Detection 

Number 

D2 

Detection 

Number 

D2 

1 

8 

1.2 

7 

msm 

» 

2 

8 

3.1 

7 

Igl 

9 

Wmm 

3 

7 

6.3 

— 

m 

— 

B 

An  alternate  strategy  is  to  always  pair  a  detection  with  a  track  if  there  is  only  one  correlation  with 
a  track.  As  before,  ambiguities  are  removed  by  using  the  smallest  statistical  distance.  Thus  track  3  in 
the  example  is  updated  by  detection  7,  track  1  is  updated  by  detection  8,  and  track  2  is  updated  by 
detection  9. 

Radar  Netting 

There  are  many  ways  of  integrating  radar  data  from  multiple  radars  into  a  single  system  track  file. 
The  type  of  radar  integration  that  should  be  used  is  a  function  of  the  radar's  performance,  the  environ¬ 
ment,  and  whether  the  radars  are  colocated  or  not.  Several  integration  methods  are: 

L  Track  Selection— Generate  a  track  with  each  radar  and  choose  one  of  the  tracks  as  the  system 

track. 

2.  Average  Track— Generate  a  track  with  each  radar,  and  weight  the  individual  tracks  to  form  a 
system  track. 

3.  Augmented  Track— Generate  a  track  with  each  radar,  choose  one  of  the  tracks  as  the  system 
track,  and  use  selected  detections  from  the  other  radars  to  update  the  system  track. 

4.  Average  Detection— Tot  a  given  time  interval  average  all  radar  detections  on  single  target  and 
use  average  detection  to  update  the  system  track. 

5.  Detection-to-  Track— Use  all  radar  detections  to  update  the  system  track;  tracks  may  or  may  not 
be  initiated  using  all  detections  from  all  radars. 

Theoretically,  the  detection-to-track  method  of  integration  yields  the  best  tracks  because  no  infor¬ 
mation  is  lost.  However,  the  detections  must  be  weighted  properly  and  care  must  be  taken  to  prevent 
bad  data  from  corrupting  good  data. 

The  fundamental  work  on  colocated  radar  integration  has  been  performed  for  the  U.S.  Navy 
where  a  typical  naval  ship  has  both  3-D  and  2-0  radars  within  several  hundred  feet  of  one  another.. 
Various  radar  integration  techniques  have  been  investigated;  however,  the  one  favored  by  the  U.S. 
Navy  is  the  detection-to-track  integration  philosophy.  Presently,  the  U.S.  Navy  is  placing  an  Integrated 
Automatic  Detection  and  Track  (IADT)  system  (called  SYS-1  120])  aboard  the  DDG-2/15  class  of 
ships  and  is  putting  a  similar  IADT  system  (SYS-2)  aboard  the  New  Threat  Upgrades  (NTU)  ships. 

When  exchanging  tracking  information  between  noncolocated  sites,  the  methods  used  are  a  func¬ 
tion  of  whether  the  sites  are  fixed  or  mobile  and  what  communication  links  are  available.  There  are 
numerous  programs  for  netting  sea,  land,  and  air  based  radars:  e.g.,  Naval  Tactical  Data  Systems 
(NTDS),  ROCC  (joint  U.S.-Canadian  system,  an  upgrade  of  SAGE  system).  Netted  Radar  Program 
(NRP),  and  various  programs  involving  airborne  sensors  such  as  AW  ACS,  NIMROD,  and  the  E2-C. 
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The  two  major  problems  associated  with  multisite  operation  are  1.  how  can  one  automatically  obtain 
gridlock  (i.e.,  properly  locating  and  aligning  radars)  using  targets  of  opportunity  and  2.  what  should  one 
transmit  over  the  communication  link  given  that  the  available  bandwidth  is  not  sufficient  to  transmit  all 
the  tracking  data. 

CONCLUDING  REMARKS 

Though  there  are  universal  tracking  systems  that  can  work  well  with  a  variety  of  radar  systems, 
there  is  no  such  thing  as  a  universal  detector  that  will  work  well  with  a  variety  of  radars  in  different 
environments.  Each  detector  should  be  designed  specifically  for  the  radar  and  the  environment  in 
which  it  is  to  operate.  That  is,  the  automatic  detector  should  be  considered  an  integral  part  of  the  radar 
system. 
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